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SPECTRAL GAMMA RAY LOGGING- WHILE-DRILLING SYSTEM 

This invention is directed toward the measure of gamma radiation from earth 
formation penetrated by a well borehole. More particularly, the invention is directed 
5 toward the determination of concentrations of naturally occurring radioactive elements in 
earth formation by analysis of energy spectra measured by at least one gamma ray 
detector while the borehole is being drilled. 

BACKGROUND OF THE INVENTION 

10 

The measure of naturally occurring gamma radiation as a function of depth within 
a well borehole is the basis of one of the earliest borehole geophysical exploration 
system. This type system, commonly referred to as a natural gamma ray logging system, 
typically comprises at least one gamma ray detector housed in a downhole tool that is 

15 conveyed along the borehole. 

One type of natural gamma ray logging system comprises a logging tool that is 
responsive to total gamma radiation emitted by the earth formation, and the tool is 
conveyed along the borehole by means of a wireline. This "total" natural gamma ray 
wireline logging system was the first type of gamma ray measurement used in borehole 

20 geophysical exploration. Since most shales are relatively rich in naturally occurring 
radioactive elements, these logs are used primarily to delineate shale from other 
formations, or used to measure the shale content of formations. This wireline logging 
system is used only after the borehole has been drilled. 

It is often advantageous to measure total natural gamma radiation while the 

25 borehole is being drilled. This is accomplished by conveying the tool along the borehole 
by means of a drill string. This type of system is commonly referred to as a total natural 
gamma ray logging-while-drilling (LWD) system. 

Yet another type of natural gamma ray logging system comprises a logging tool 
that measures a spectrum of gamma radiation emitted by the earth formation. The 

30 spectrum is defined a measure of intensity of radiation as a function of radiation energy. 
This type of logging system is commonly referred to as a spectral gamma ray logging 
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system. Spectral gamma ray logging tools are typically conveyed along the borehole by 
means of a wireline. Low count rate and detector stabilization are major problems in any 
type of natural spectral LWD systems 

Most naturally occurring gamma radiation found in earth formations is emitted by 
5 potassium (K) and elements within the decay chains of uranium (U) and thorium (Th). 
Energy of naturally occurring gamma radiation measurable in a borehole environment 
typically spans a range of about 0.1 to less than 3.0 million electron Volts (MeV). The 
elements K, U and Th emit gamma radiation at different characteristic energies. 
Components of radiation from K, U and Th contributing to the total measured gamma 

10 radiation can, therefore, be obtained by identifying these characteristic energies using 
spectral gamma ray logging system. Through system calibration and modeling, these 
components can be subsequently related to the corresponding elemental concentrations of 
these elements within the formation. Elemental concentrations of K, U and Th can be 
used to determine parameters in addition to shale content obtained from total natural 

15 gamma ray logs. These additional parameters include, but are not limited to, clay typing, 
lithology identification, fracture detection, and radioactive tracer placement. 

As in all nuclear logging systems, statistical precision of a measurement is 
maximized when the count rate of the radiation detector used to obtain the measurement 
is maximized. Naturally occurring gamma radiation is typically much less intense than 

20 gamma radiation induced in formation materials by sources of radiation within a logging 
tool. It is important, therefore, to design natural gamma ray logging tools to maximized 
measured gamma radiation count rate. 

Measured count rate can be optimized by designing tool housings (both total 
gamma ray and spectral gamma ray) so that gamma radiation attenuation within the 

25 housing is minimized. The lower energy region of the measured spectrum is especially 
important in spectral gamma ray logging systems. Wireline spectral gamma ray logging 
tools often employ a tool housing fabricated with material of relatively low atomic 
number, rather than heavier (and stronger) materials such as steel. These so called "low 
Z" tool cases minimize gamma ray attenuation, especially at the lower end of the energy 

30 spectrum, thereby maximizing measured count rate for a given radiation intensity and 
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detector size. Low Z materials often do not meet structural requirements of LWD 
systems. 

Measured count rates can further be maximized through tool detector design. Due 
to the relatively high energies of the characteristic K, U, and Th gamma radiation, it is 
5 advantageous for the gamma ray detector of a given type to be dimensioned as large as 
practically possible to react with, and thereby respond to, these radiations. Typically, 
larger detectors can be disposed in wireline tools with less attenuating material between 
the detector and the formation. LWD systems employ a relatively thick tool housing, 
which is typically a collar with a drilling fluid flow conduit passing through the collar. 

10 A gamma ray detector comprising a scintillation crystal and a cooperating light 

sensing device, such as a photomultiplier tube, typically yields the highest spectral 
gamma ray detector efficiency for a given detector volume. Gamma ray detectors 
undergo significant temperature changes during a logging operation. The gain of a 
photomultiplier tube changes as the temperature and, to a lesser extent, counting rate 

15 changes. Gain changes, often referred to as gain "shifts", adversely affect gamma ray 
spectral analysis. Typically, a 100 degree Centigrade (°C) change in temperature causes 
100% change in gain. Temperature variations of this order of magnitude are not 
uncommon in wireline or LWD logging operations. It is, therefore, necessary to 
compensate for detector gain changes in order to obtain accurate and precise spectral 

20 gamma ray measurements. This compensation is especially difficult to achieve in LWD 
systems. As an example, significant gain changes can occur over a relatively short time 
interval. The data rates of available LWD telemetry systems between the downhole tool 
and surface equipment are typically too low to effectively monitor and to correct for 
rapidly occurring gain shifts. Automatic downhole gain control is, therefore, highly 

25 desirable in LWD systems. 

As mentioned previously, naturally occurring gamma ray spectral measurements 
are typically low count rate. It is, therefore, desirable to use as much of the measurable 
gamma ray spectrum as possible in order to maximize statistical precision. Shock and 
vibration effects on low count rate systems can distort spectral shape, especially at the 

30 lower energy region of the measured spectrum. This problem is especially prevalent in 
LWD systems, which are exposed to harsh drilling environments. 
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SUMMARY OF THE INVENTION 

This present invention is directed toward a spectral gamma ray logging-while- 
5 drilling (LWD) system. The system is designed to yield elemental concentrations of 
naturally occurring radioactive material such as K, U and Th. It should be understood, 
however, that the system can be used to obtain spectral measurements of any type of 
gamma radiation encountered in a borehole environment. 

The LWD downhole assembly or "tool" comprises a drill collar that is attached to 

10 the lower end of a drill string. A drill bit terminates the lower end of the tool. Sensor, 
electronics and downhole telemetry elements are disposed within the collar. The tool is 
conveyed along a well borehole by means of the drill string, which is operated by a rotary 
drilling rig at the surface of the earth. Information from the tool is telemetered to the 
surface via a telemetry link and received by a surface telemetry element contained in 

15 surface equipment that is operationally attached to the drilling rig. Information can also 
be transmitted from the surface equipment to the tool via the telemetry link. 

The sensor element comprises one or more gamma ray detectors that are disposed 
as close as practical to the periphery of the tool. This minimizes intervening material 
between the one or more detectors and the source of gamma radiation, which is earth 

20 formation penetrated by the borehole. As a result of this detector geometry, spectral 
degradation is minimized and measured count rate is maximized for a given detector size. 
The detector geometry also allows an azimuthal spectral gamma ray measurement in a 
plane essentially perpendicular to the axis of the tool. The one or more gamma ray 
detectors preferably comprises a scintillation crystal optically coupled to a light sensitive 

25 device such as a photomultiplier tube. The detector element is calibrated under known 
conditions and at a "standard" detector gain. The sensor element can also contain a 
system, such as a magnetometer, that senses the orientation of the tool within the 
borehole. 

Output signals from the sensor element are input to the electronics element. The 
30 signals are amplified using appropriate preamplification and amplification circuits. 
Amplified sensor signals are then input to a processor for subsequent processing. High 
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voltage for the one or more gamma ray detectors is provided by an adjustable high 
voltage power supply within the electronics element. Changes in temperature or, to a 
lesser extent, changes in measured gamma ray count rate result in detector gain change. 
Peak structure location and continuum regions of measured gamma ray spectra are 
5 monitored by the processor. Any gain change is detected using predetermined 
relationships and criteria stored within the processor. A gain correction signal 
representative of the magnitude of the gain change is generated by the processor and 
input to the adjustable high voltage power supply thereby adjusting detector high voltage 
such that the gain is restored to the standard gain. This gain control system is automatic, 

10 and requires no intervention from the surface. 

With detector gain stabilized to standard gain, elemental concentrations of K, U 
and Th are determined in the processor using predetermined relationships. These 
elemental concentrations can be input to the downhole telemetry element and telemetered 
to the surface. Alternately, gain stabilized spectral data can be input to the downhole 

15 telemetry element and telemetered to the surface for subsequent processing. Spectral 
gamma ray data and elemental concentration determinations can be recorded by a data 
storage means within the electronics element, and subsequently extracted for processing 
and analysis when the tool is returned to the surface of the earth. 

Elemental concentrations of K, U and Th are determined as a function of depth as 

20 the tool is conveyed along the borehole. If a plurality of gamma ray detectors is used, the 
gain adjusted spectral responses of the detectors are combined to obtain the desired 
elemental concentrations. Preferably the detector responses are combined prior to 
computation of elemental concentrations. 

The peripheral detector geometry also allows an azimuthal spectral gamma ray 

25 measurement and corresponding azimuthal elemental concentration determinations in a 
plane that is essentially perpendicular to the axis of the tool. Azimuthal reference is 
obtained by using a tool orientation sensitive device such as a magnetometer disposed 
within the sensor or electronics element. If a single detector is used, azimuthal 
measurements can be obtained only when the tool is being rotated by the drill string. A 

30 plurality of detectors yields azimuthal information when the tool is rotating or "sliding" 
along the borehole without rotating. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

So that the manner in which the above recited features, advantages and objects the 
5 present invention are obtained and can be understood in detail, more particular 
description of the invention, briefly summarized above, may be had by reference to the 
embodiments thereof which are illustrated in the appended drawings. 

Fig. 1 illustrates the spectral natural gamma ray LWD system displayed as a 

whole; 

10 Fig. 2 is a functional diagram of major elements and components of the spectral 

LWD system; 

Fig. 3a is a cross sectional view of a spectral LWD tool sensor element 
comprising one gamma ray detector; 

Fig. 3b is a side sectional view of the sensor element comprising one gamma ray 
15 detector; 

Fig. 4 is a cross sectional view of a spectral LWD tool sensor element comprising 
three gamma ray detectors; 

Fig. 5 is a typical gamma ray natural gamma ray spectrum measured with the 
spectral gamma ray LWD tool; 
20 Fig. 6 shows a relationship between the slope of the Compton region of the 

measured gamma ray spectrum and the temperature of the gamma ray detector; 

Fig. 7 illustrates relationships between detector gain adjustment factor required to 
obtain standard detector gain, a required high voltage adjustment to obtain standard 
detector gain, and detector temperature; 
25 Fig. 8 is a more detailed view of a measured gamma ray spectrum illustrating how 

a spectral peak position is used to obtain a second order detector gain adjustment; 

Fig. 9 is a graphical illustration of a method for locating statistically significant 
peak structure in a measured gamma ray spectrum; 

Fig. 10 is a flow chart showing steps for automatically controlling the gain of a 
30 gamma ray detector using a measured spectral analysis method; 
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Fig. 1 1 illustrates basic concepts used for automatically controlling the gain of a 
gamma ray detector using a detector source gain control method; 

Fig. 12 is a flow chart showing steps for automatically controlling the gain of a 
gamma ray detector using the detector source gain control method; 
5 Fig. 13 shows additional hardware components needed to obtain azimuthal natural 

gamma ray measurements using the LWD system; 

Fig. 14 is an example of a spectral gamma ray LWD log showing concentrations 
of K, U and Th as a function of depth within a borehole; and 

Fig. 15 is an example of an azimuthal spectral gamma ray LWD log showing 
10 concentrations of K, U and Th as a function of azimuth around the borehole and as a 
function of depth within the borehole. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

15 Details of the preferred embodiments of the LWD spectral gamma ray logging 

system are presented in sections. System hardware is first disclosed. This is followed by 
disclosure of methodology used to monitor measured gamma ray spectra, and to stabilize 
the gain of these spectra as borehole temperature varies. Two gain stabilization methods 
are disclosed. With both, stabilization is accomplished in real time and without operator 

20 intervention. Once gain stabilization has been obtained, methods for determining 
elemental concentrations of naturally occurring K, U and Th are discussed. Finally, 
measures of total and azimuthal concentrations of K, U and Th are discussed, and "log" 
presentations of these measurements are illustrated. 

The invention is directed toward the measure of gamma radiation that occurs 

25 naturally in earth formation. It should be understood, however, that the basic concepts of 
the invention are applicable for quantitative measurements of any type of gamma 
radiation wherein one or more gamma ray detectors are subject to gain shifts. 

HARDWARE 

30 Fig. 1 illustrates the LWD system 15 displayed as a whole. A downhole assembly 

or "tool" comprises a drill collar 10 that is attached to the lower end of a drill string 18. 
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A drill bit 1 1 terminates the lower end of the collar 10. Within the collar 10 are disposed 
a sensor element 12, an electronics element 14, and a downhole telemetry 16. The tool is 
conveyed along a well borehole 20, defined by borehole walls 21 and penetrating 
formation 22, by means of the drill string 18. The drill string 18 is operated from the 
5 surface 24 of the earth by rotary drilling rig, which is only illustrated conceptually at 26 
since such rigs are well known in the art. 

Information from the tool is telemetered to the surface of the earth 24 via a 
telemetry link (illustrated conceptually by the arrow 23) and received by a surface 
telemetry element (not shown) contained in surface equipment 28 that is operationally 

10 connected to the drilling rig 26. Information can also be transmitted from the surface 
equipment 28 to the tool via the telemetry link 23. 

More details of the sensor element 12, the electronics element 14. and the 
downhole telemetry element 16 and their operating relationships are shown in the 
functional diagram of Fig. 2. The sensor element 12, illustrated conceptually as a broken 

15 line box, comprises at least one gamma ray detector comprising a scintillation crystal 30 
and an optically coupled photomultiplier tube 32. Output signals from the 
photomultiplier tube are input to the electronics element, whose components are enclosed 
by the broken line box designated as 14. The signals are amplified using appropriate 
preamplification and amplification circuits 34. Amplified sensor signals are input to a 

20 processor 38. Voltage for the photomultiplier tube 32 is provided by an adjustable high 
voltage power supply 36 within the electronics element 14. 

Still referring to Fig. 2, the processor 38 provides means for automatically 
controlling the gain of the at least one gamma ray detector, and is also used to process 
signals from the gamma ray detector to obtain elemental concentrations of K, U and Th. 

25 As mentioned previously, the tool is preferably calibrated at the surface to a "standard" 
gain. While logging, the temperature of the tool and elements therein change. Changes 
in temperature or, to a lesser extent, changes in measured gamma ray count rate result in 
detector gain change. Gain changes are reflected in the energy spectrum of the measured 
detector signals. Measured features of the spectrum are used to correct for these gain 

30 changes, as will be discussed in detail in subsequent sections of this disclosure. A gain 
correction signal representative of the magnitude of the gain change is generated by the 
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processor 38 and input to the adjustable high voltage power supply 36 thereby adjusting 
detector high voltage so that the gain is restored to the "standard" gain. Elemental 
concentrations of K, U and Th are determined from gain corrected detector spectra in the 
processor 38 using predetermined relations as will be discussed subsequently. Elemental 
5 concentrations of K, U and Th are input to the downhole telemetry element 16 and 
telemetered, via the telemetry link 23, to the surface telemetry element contained in 
surface equipment 28. 

Figs. 3a and 3b are cross sectional and side sectional views, respective, of the 
collar 10 in the region of the sensor element and depict a sensor element comprising one 

10 gamma ray detector. Fig. 3a shows the cross section A- A' of the collar 10 with the axis 
of a drilling flow conduit 44 displaced from the axis of the collar. A detector channel on 
the periphery of the collar 10 and defined by the surfaces 40 receives the gamma ray 
detector comprising a scintillation crystal 30 such as Nal, Csl, BGO and the like. The 
scintillation crystal 30 is encapsulates in a hermetically sealed, light reflecting casing 42. 

15 The volume 46 is preferably filled with a material such RTV, epoxy and the like. The 
collar 10 is surrounded by a thin sleeve 48 in the region of the sensor element 12. Fig 3b 
shows the side section B-B' which includes the major axis of the collar 10. A 
photomultiplier tube 32 is optically coupled to the scintillation crystal 30. Electrical 
leads to the photomultiplier tube are not shown for purposes of clarity. 

20 Again referring to both Figs 3a and 3b, it is apparent that the scintillation crystal 

30 is disposed as close as practical to the periphery of the collar 10. This minimizes 
intervening material between the detector and the source of gamma radiation, which is 
earth formation penetrated by the borehole (not shown). By displacing the axis of the 
flow conduit 44, the diameter, and thus the efficiency, of the detector is maximized. For 

25 typical LWD equipment, the diameter of the scintillation crystal can be 2 inches (5.1 
centimeters) or larger and still maintain structural specifications of the collar 10. As a 
result of this detector geometry, gamma ray spectral degradation is minimized and 
measured count rate is maximized for a given detector size. 

Fig. 4 illustrates a sensor element comprising three gamma ray detectors. This 

30 cross section view shows the scintillation crystals 30 of each detector. Each crystal 30 is 
encapsulated in hermetically sealed, light reflecting casing 42, and is disposed in a 
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detector channel defined by the surfaces 40. The channels are arranged at 120 degree 
angular spacings. The collar 10 is in the region of the sensor element and is again 
surrounded by a thin sleeve 48. A side section view has been omitted for brevity, but a 
photomultiplier tube (not shown) is again optically coupled to each scintillation crystal 
5 30. As with the single detector sensor element shown in Figs 3a and 3b, it is apparent 
that the scintillation crystals 30 are disposed as close as practical to the periphery of the 
tool thereby minimizing intervening material between the detectors and the source of 
gamma radiation within the earth formation. Using the multiple detector configuration, 
the axis of the flow conduit 44 is. coincident with the axis of the collar 10. For typical 

10 LWD equipment, the diameter of each scintillation crystal is limited to about 1.5 inches 
(3.8 centimeters) so that structural specifications of the collar can be maintained. 
Compared with the single scintillation crystal sensor element configuration shown in Figs 
3a and 3b, the efficiency of each detector on Fig. 4 is reduced. Signals from each 
detector can, however, be combined to obtain a total sensor element efficiency that equals 

15 or exceeds the efficiency of the single detector configuration. In addition, the three 
detector sensor element configuration offers advantages in azimuthal spectral gamma ray 
measurements that will be discussed in a subsequent section of this disclosure. 

It should be understood that the multiple detector sensor element configuration is 
not limited to the three detector configuration shown in Fig. 4. It should also be 

20 understood that angular spacing between the multiple detectors need not be equal. 
Finally, it should be understood that the dimensions of the multiple detectors need not be 
the same. 

GAIN STABILIZATION 

25 Two methods of detector gain stabilization are disclosed. The first method will be 

referred to as the "measured spectral analysis" method, and the second method will be 
referred to as the "detector source" gain correction method. The gain of an LWD gamma 
ray detector can changfc significantly and rapidly in the harsh borehole drilling 
environment. Telemetry links between the tool and the surface are relatively slow and do 

30 not permit gain monitoring and correction from the surface. Gain control must be 
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implemented automatically within the tool. Both of the disclosed methods can be used to 
effectively control gamma ray detector gain. 

Considering the importance of gain control and the harshness of the borehole 
environment, it is desirable to use both methods. The two methods can be used with one 
5 serving as a primary method for gain control, and the second serving as a back-up method 
for gain control. Alternately, both methods can be used simultaneously, and the results 
combined to obtain a detector gain correction. Such a combination can take a variety of 
forms including a simple numerical average or a weighted average. 

10 Measured Spectral Analysis Method 

Fig. 5 is a typically natural gamma ray spectrum measured in earth formation with 
a scintillation type gamma ray detector. The spectrum comprises measured gamma ray 
intensity as a function of gamma ray energy, represented by the curve 50. The abscissa is 
gamma ray energy in million electron Volt (MeV), and the ordinate in the natural log of 

15 measured count rate per increment of energy. The increments of energy are represented 
energy channels or "channels" on the top scale abscissa. Representative peak structure 
from K, U and Th is shown at energies 1.46 MeV, 1.76 MeV and 2.61 MeV, respectively. 
During tool calibration, the detector high voltage is adjusted to give a detector gain for 
which specific energies of gamma radiation fall within predetermined energy channels. 

20 This gain is referred to as the "standard" gain. Tool calibration will be discussed further 
in a subsequent section of this disclosure. 

The Compton scatter region of the spectrum comprises formation gamma 
radiation that has undergone several collisions in intervening material before it reaches 
the gamma ray detector. This region of the spectrum is identified by the numeral 52, and 

25 terminates at the low energy region by the "hump" 54 at energy 56. This exponential- 
looking region 52 contains no direct contributions from K, U, and Th gamma radiation. 
The slope of this region is mainly a function of the photomultiplier tube gain, and can be 
used to monitor detector gain. 

Fig. 6 is a plot of the measures of the slope of the Compton region 52 as a 

30 function of detector temperature in degrees Centigrade (°C). Measured values of slope 
(ordinate) versus temperature (abscissa) are indicated by the data points 62. A curve 60 

11 



AES 04-001 



fitted through the data points 62 shows that there is essentially 1 : 1 correlation between 
the slope and detector temperature. 

Curve 66 of Fig. 7 shows the relationship between a multiplicative first order gain 
adjustment factor Fj and corresponding detector temperature, where 

5 

(1) G st d = FiG 0 bs . 

Gobs is the observed detector gain, and G st d is the previously defined "standard" gain for 
which the tool is calibrated. Relationships of slope versus temperature shown graphically 

10 as curve 60 on Fig. 6 can be combined with gain adjustment factor as a function of 
temperature shown graphically as curve 66 in Fig. 7 to eliminate the temperature 
parameter. This combination yields a functional relationship between the measured slope 
and the desired first order gain adjustment factor Fi. 

At this point, detector high voltage can be adjusted to correct detector gain for 

15 temperature effects. Curve 68 of Fig 7 shows the relationship between required high 
voltage adjustment to obtain G st d and detector temperature. Once Fi has been obtained as 
described above, the high voltage Vi required to obtain G st d can then be determined. The 
following example is presented as a graphical solution. Assume that from a measure of 
the Compton slope, it has been determined that Fi = 1.1 as indicated in Fig. 7 at 70. A 

20 horizontal line is projected until it intersects the curve 66 at a point 71. A vertical line is 
projected until it intersects the curve 68 at point 72. Finally, a horizontal line is projected 
to the right ordinate at point 73 giving a required high voltage correction of +10 volts. 
Referring again to Fig. 2, 10 volts added to the photomultiplier tube 32 by the adjustable 
high voltage power supply will adjust the gain of the gamma ray detector to the standard 

25 gain G std . 

To summarize, the relationships shown graphically in Fig. 6 and Fig. 7 are 
combined to develop a relationship between measured slope in the Compton region as a 
function of high voltage required to maintain standard detector gain. The example 
discussed above is graphical, but it should be understood that the solution can be reduced 
30 to analytical form suitable for computation in a processor. A measure of slope can, 
therefore, be used to correct for detector gain shifts. Using this slope to predict gain 
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changes is extremely robust since it is calculated from a number of data points, it is 
immune to electronics noise, and it has been found that is not highly affected by borehole 
and formation conditions. A first order detector gain correction can be made for the most 
severe gain changes ranging from about -60% of the standard gain to +150% of the 
5 standard gain. Such severe gain changes are induced by equally severe changes in 
temperature ranging from about -60 °C to about +150 °C, with 25 °C being the "standard" 
calibration temperature. 

As discussed above, a measure of slope can be used to obtain a voltage 
adjustment required to obtain standard gain. This is a first order correction considering 

10 the magnitude of gain changes that can be handled. Second and third order corrections 
are made to further increase precision and accuracy of the detector gain settings. Because 
of these additional corrections, it is preferred not to adjust detector voltage after the first 
order gain correction. Instead, the spectrum is adjusted by adjusting the count rate per 
energy channel using the first order gain correction factor Fi. Using the previous 

15 example of Fi = 1.1, the width of each channel is "stretched" by 10 percent, and the 
measured count rates are redistributed over the wider channels. This methodology can be 
thought of a "software" gain shift, and is preferably performed in processor 38. The 
detector gain adjustment is now further refined by examining a predominate peak in this 
modified energy spectrum. Fig. 8 shows a portion of the full spectrum 50 shown in Fig. 

20 5, and includes a peak 80 at 1.46 MeV from K. This peak is typically the most prominent 
peak in the spectrum, and is suitable for the second order gain adjustment. After 
modifying the spectrum using methodology previously discussed, it is observed that the 
maxima of the peak 80 falls in a channel Chobs at 78. Tool calibration requires that the 
energy corresponding to this peak maxima fall in a modified channel Ch st d at 76, where 

25 the channels have been adjusted in width by the first order correction Fj. The gain of the 
detector is further adjusted to a second order correction so that the maxima of the peak 80 
falls in Ch s td. This is again accomplished with a "software" gain shift by adjusting the 
widths of the energy channels and redistributing the measured count rates to form a 
second modified spectrum. A second order gain correction F 2 is obtained from the 

30 equation 



13 



AES 04-001 



(2a) F 2 = (Chstd/Chobs) . 

The gain correction process can be terminated at this point, with no further refinement of 
the gain setting. If this option is chosen, the corresponding high voltage setting V 2 
5 required to obtain this second order corrected standard gain and is expressed 
mathematically as 

(2b) V 2 = F 2 V, 

10 Peaks in measured gamma ray spectra such as the curve 50 shown in Fig. 5 are 

identified with respect to channel (and corresponding energy) using a second derivative 
algorithm after a strong filter is applied. Locations of the peaks are determined using a 
Guassian fit around each peak. Fig. 9 illustrates the effectiveness of the detection and 
location method. Curve 86 represents the second derivative d 2 C/dCh 2 (right ordinate) of 

15 spectral count rate C as a function of corresponding energy channel Ch for a spectrum of 
the type shown in Fig. 5. Peaks are indicated when the curve 86 crosses d 2 C/dCh 2 = 0. 
Since the spectrum exhibits statistical variations, the curve can also exhibit a "zero 
crossing" due to statistics of the count rate rather than actual peak structure. It is, 
therefore, necessary to identify a non-statistical or "noise" zero crossings from a true zero 

20 crossing which indicate a peak. The curves 82 and 84 represent upper and lower limits of 
standard deviation of the count rate C as a function of energy channel, and are presented 
in arbitrary count rate units of the left ordinate. Only zero crossing excursions that 
extend outside the standard deviation "envelope" are considered as statistically significant 
indications of a peak. There is only one such indication in the curve 86 at a channel 

25 identified by the numeral 88. This corresponds to the K peak at 1.46 MeV shown clearly 
in Fig. 5. 

It is preferred to still further refine the gain setting. Once the second modified 
energy spectrum has been computed, all statistically significant peaks in the second 
modified spectrum are located using the peak location technique discussed above. Once 
30 the energy channels in which these peaks are observed, they are compared with their 
corresponding "standard" channels determined in tool calibration. This methodology is 
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similar to the single-peak methodology used for the second order gain correction, but all 
peaks are used in this third order correction. Channel widths are again adjusted and count 
rates redistributed so that all identifiable peaks fall in their corresponding standard energy 
channels, corrected for the first and second order software gain adjustments. This is the 
5 third order gain correction and yields a third order gain adjustment factor F3. The 
uncorrected detector voltage V is now adjusted to obtain a corrected voltage V cor using 
the relationship 

(3a) V cor = F 1 F 2 F 3 V . 

10 

The corrected detector gain G cor is 
(3b) G cor = HV cor 

15 where H is a multiplicative constant relating the third order corrected high voltage V cor to 
the fully corrected detector gain G cor . Channel widths are reset to their original values. 

The measured spectral analysis method of automatic gain correction is 
summarized in the flow chart shown in Fig. 10. The slope of the Compton region 52 (see 
Fig. 5) is measured at 90. The first order gain correction Fi is determined at 92 using the 

20 measured slope as discussed above. An identifiable peak is located in the measured 
spectrum at 94. The second order gain correction F2 is determined at 96 using previously 
discussed methodology and equation (2a). All statistically significant peaks in the 
measured gamma ray spectrum are determined at 98. The third order gain correction 
factor F3 is obtained at 100 by a software gain adjustment that positions all peaks in 

25 alignment with their assigned energy channels obtained at tool calibration. The 
corresponding high voltage V cor required for standard detector gain is also determined at 
100 using equation (3a). The corrected detector gain G cor is set at 102 using V 3 and 
equation (3b). It again noted that an actual voltage adjustment is made only after the 
third order correction, with software gain adjustments being used in the first and second 

30 order corrections. It should also be understood that other algorithms can be used to 
obtain suitable software gain adjustments for the first and second order corrections. 
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Detector Source Gain Correction Method 

Gamma ray detector gain can be monitored using an alternate technique. A small 
radioactive "detector" source is disposed near or within the one or more scintillation 
5 crystals comprising the natural gamma ray LWD sensor element. The detector source 
generates a "calibration" peak in the measured gamma ray spectrum. If the gain of the 
measured spectrum changes, the position of the calibration peak shifts with the change in 
gain. A measure of position of the calibration peak can, therefore, be used to monitor and 
to correct detector gain. 

10 The calibration peak is preferably at a relatively low energy so that it will not 

interfere with higher energy radiation from K, U and Th used to determine elemental 
concentrations. A suitable detector source is Americium-241 ( 24I Am) which emits 
gamma radiation at 0.060 MeV. Referring again to Fig. 5, a typical measured natural 
gamma ray spectrum encompasses an energy range of 0.0 to 3.0 MeV over typically 256 

15 channels. With this "standard" detector gain setting, the low energy end of the spectrum, 
which includes the 0.060 MeV calibration peak, is very susceptible to electronics noise. 
In addition, since the spectrum spreads over about 3.0 MeV and typically 256 energy 
channels, the 241 Am peak occupies only about 3 out the 256 channels, which makes it 
difficult to precisely locate the peak position. 

20 Attention is directed to Fig. 1 1, which conceptually depicts the low energy region 

of a gamma ray spectrum measured with a detector comprising a 241 Am detector source. 
The curve 1 10 shows the peak structure with the detector set at standard gain. Since the 
peak occupies only three energy channels, it is difficult to locate the center of the peak 
using previously discussed methods. This peak can, in principle, be used as shown to 

25 stabilize detector gain. Any gain stabilization using the 241 Am peak at 0.060 MeV would, 
therefore, be subject to large errors, especially at the higher energy region of the spectrum 
used for elemental concentration calculations. 

The precision of detector gain stabilization using a low energy detector source and 
calibration peak is improved using dual gain circuitry. The spectrum signal from the 

30 gamma ray detector is "branched" and input into first and second amplification circuits 
comprising the dual gain circuitry. The first circuit comprises a standard amplification 

16 



AES 04-001 

circuit and generates a spectrum with standard gain. The 0.060 MeV 241 Am peak in this 
standard spectrum is shown at 110 in Fig. 11. For purposes of gain stabilization, the 
spectrum signal is input to a second amplification circuit with a gain factor of N greater 
than that of the standard gain. The second amplification circuit generates an amplified 
5 spectrum with an "amplified" gain. For purposes of discussion, it will be assumed that N 
= 10, although it should be understood that other values of N can be used. Curve 1 12 is 
the amplified gain spectrum showing the 241 Am peak at 0.060 MeV amplified by a factor 
of N = 10. The peak now occupies about 30 energy channels, and previously discussed 
peak location methods are used to determine that the center of the peak is in energy 
10 channel P 0 bs, which identified at 114. From tool calibration, it is known that energy 
0.060 MeV should fall in energy channel P st d for the standard gain, or in channel N x P st d 
for the amplified gain, shown at 1 16 in Fig. 11. The detector high voltage V is adjusted 
to a corrected value, V cor , using the relationship 

15 (4) V cor = V(NP $ A s ) . 

A signal proportional to (N P s td/P 0 bs) is preferably generated in the processor 38 and input 
to the adjustable high voltage power supply 36. This generates the corrected high voltage 
V cor supplied to the detector. Corrected standard gain amplification G cor is expressed by 
20 the relationship 

(5) Gcor = H V cor 

where H, as in equation (3), is a multiplicative constant relating high voltage to detector 
25 gain. 

It should be understood that various methods can be used to increase the detector 
gain by a factor of N. As an example, the amplification circuit 34 (see Fig. 2) can contain 
a high gain element which, under the control of the processor 38, "branches" detector 
input. The gain of one branch of the input by a factor on N thereby forming the amplified 
30 gain of the calibration peak for purposes of calibration source gain stabilization. 
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The gain stabilization method using a radioactive "detector" source is summarized 
in the flow chart of Fig. 12. The low energy portion of the measured gamma ray 
spectrum is increased by a factor of N at 120. The energy channel P 0 bs, in which the 
stabilization peak is maximum, is determined at 122 using a suitable peak location 
5 technique. High voltage V cor required to position the peak maxima in energy channel N x 
P st d is determined (see equation (4)) at 124, and the correct detector gain G cor is set at 126 
(see equation (5)). 

Hybrid Correction Method 

10 The measured spectral analysis method and the detector source gain correction 

methods can be combined to yield a hybrid gain control method. A calibration source is 
disposed in within or in the immediate vicinity of at least one gamma ray detector. When 
in the borehole, this detector produces a gamma ray spectrum comprising a first 
component from naturally occurring radioactive elements within the formation and a 

15 second component from the calibration source. A first detector gain correction is 
determined from spectral features from the naturally occurring gamma radiation as 
previously discussed in the measured spectral analysis method. A second gain correction 
is determined from the calibration source component as previously discussed in the 
detector source gain correction method. The first and second gain corrections are 

20 combined to correct for gain shift of the detector. 

ELEMENTAL CONCENTRATION DETERMINATIONS 

With detector gain stabilized to "standard" gain, elemental concentrations of K, U 
25 and Th are determined, preferably in the processor 38 of the electronics element 14 (see 
Figs. 1 and 2), from measured spectral data. These elemental concentrations can be input 
to the downhole telemetry element 16 and telemetered via the telemetry link 23 to the 
surface equipment 28. Alternately, the spectral data can be input to the downhole 
telemetry element 16 and telemetered to the surface equipment 28 for subsequent 
30 processing. Since the telemetry bandwidth is limited and the gamma ray spectra are 
much more data intensive than the elemental concentrations determined therefrom, it is 
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preferred to telemeter the elemental concentrations of K, U and Th to the surface. 
Alternately, spectral gamma ray data and elemental concentration determinations can be 
recorded by a data storage means within the electronics element, and subsequently 
extracted from processing and analysis when the tool is returned to the surface. 
5 The following methodology is preferred for determining elemental concentrations 

of K, U and Th. It should be understood, however, that other spectral processing 
methods such as spectrum stripping, peak area analysis and the like can be used to 
determine concentrations of K, U and Th. The required elemental concentration 
calibration constants are obtained at tool calibration. 
10 Elemental concentrations are obtained by solving the matrix equation 

(6) [C] = [A][M] . 

[C] is a m x 1 column matrix comprising elements Ci (i = 1, m) representing 
15 count rate recorded in energy channel i (see Fig. 5). Typically 256 energy channels (m = 
256) are used, although more or fewer channels can be used within the scope of the 
invention. 

[A] is a m x j matrix comprising elements a^ with (i = 1, m) and (j = K, U, Th). 
Physically, the element is the sensitivity of energy channel i to the element j, typically 

20 in units of counts per second per part per million (U and Th) or counts per second per 
percent (K). The matrix [A] comprises calibration constants, is referred to as a 
"sensitivity" matrix, and is determined at tool calibration. At tool calibration, the 
response of the tool is measured in formations containing known concentrations of K, U 
and Th, and in "standard" borehole conditions, and with the one or more detectors in the 

25 sensor element operated at "standard" gain G cor . 

[M] is a j x 1 column matrix comprising elements Mj (j = K, U, Th) which are the 
parameters of interest, namely the formation elemental concentrations of K, U and Th. 
Mk is in percent, and My and Mjh in parts per million (ppm). The desired elemental 
concentrations are obtained by solving equation (6) for [M], preferably using a weighted 

30 least squares fit. 
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Measured gamma ray spectra from one or more gamma ray detectors in the sensor 
section are tracked as a function of depth of the tool in the borehole 20 (see Fig. 1). If the 
sensor element 12 comprises only one gamma ray detector as shown in Figs. 3a and 3b, 
the elements of the matrix [C] are obtained from that detector. If the sensor element 
5 comprised a plurality q of detectors, such as the q = 3 embodiment shown in Fig. 4, the 
elements of the [C] matrix are obtained by- combining responses of the q detectors, 
typically by simply summing the responses if all detectors exhibit equal sensitivity. 

AZIMUTH AL ELEMENTAL CONCENTRATION DETERMINATIONS 

10 

The spectral gamma ray LWD system can be used to measure elemental 
concentrations Mk, Mu and M-m as a function of azimuth within the borehole as well as a 
function of depth within the borehole. Azimuthal measurements require additional 
components disposed preferably within the electronics element 14. Fig. 13 is a functional 

15 diagram of components added to the electronics element shown in Fig. 2 so that 
azimuthal elemental concentrations M K , Mu and M™ can be determined. A device 
sensing tool orientation, such as a magnetometer 130, and a clock 132 are operationally 
connected to the processor 38. As in the previous discussion of Fig. 2, signals from the 
one or more gamma ray detectors and amplifier circuits are input to the processor at 136. 

20 The processor 38 again controls detector gain adjustments of the detectors at 138 through 
the adjustable high voltage power supply 36. Spectral and elemental concentrations are 
output from the processor at 134 as described below. 

As the tool rotates through 360 degrees, gamma ray spectra of the form shown in 
Fig. 5 are measured during discrete time intervals At, where these elements are defined by 

25 the clock 132 cooperating with the processor 38. Spectra are stored in bins according to 
time intervals At in which they are measured. The time interval At is preferably about 50 
milliseconds. During each time interval, the average reading of the magnetometer 130 is 
determined thereby defining an azimuth sector associated with each time interval, and 
thereby assigned an azimuth sector to each bin. Each bin, therefore, contains a gamma 

30 ray spectrum measured at a known borehole azimuth sector. The spectral data binning, 
and the averaging of magnetometer readings during each time interval, are controlled by 



20 



AES 04-001 



the processor 38. Binned spectral data and corresponding azimuth sectors are preferably 
stored in the processor 38. The process is repeated through multiple 360 degree rotations 
within a given depth interval Ad in order to maximize statistical precision of each natural 
gamma ray spectrum stored in each bin. Gain stabilization techniques previously 
5 disclosed in detail are used to control the gain of each binned spectrum. Previously 
discussed data analysis methods are used to compute the matrix [M] for each binned 
spectrum thereby yielding elemental concentrations M K , Mu and Mm for each azimuth 
sector around the borehole. 

If the sensor element comprises a plurality of detectors, detector outputs are 

10 phased by the processor 38 so that as each detector rotates through each azimuth sector, 
output from that detector is stored within the bin corresponding to that azimuth sector. 

The tool can be conveyed along the borehole without rotating. This conveyance 
is commonly referred to as "sliding". If the sensor element 12 comprises only one 
gamma ray detector, azimuthal natural gamma ray spectral measurements can not be 

15 made when the tool is sliding. If the sensor element comprises a plurality of gamma ray 
detectors, azimuthal spectral measurements can be obtained while sliding. The 
magnitudes of the azimuth sectors are determined by the number of detectors in the 
sensor element. For the sensor element comprising three detectors on 120 degree centers 
as shown in Fig, 4, each azimuth sector would be 120 degrees. This yields an azimuthal 

20 resolution that is typically inferior to that obtained with the tool rotating and with time 
intervals At of about 50 milliseconds. 

LOG PRESENTATIONS 

25 Fig. 14 shows an example 140 of a natural gamma ray LWD log presentation of 

M K , Mu and as a function of depth in the borehole. The quantities M K , Mu and Mjh 
are computed from measured spectral data [C] using equation (6). Units for 
concentrations of K (%), U (ppm) and Th (ppm) are shown in the fields 141, 143 and 
145, respectively. Scales are typically in % per chart division for K and ppm per chart 

30 division for U and Th. Concentrations of Mr, Mu and Mn, are shown as a function of 
depth 148 in the borehole by curves 142, 144 and 146, respectively. As an example, 
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excursions 147 and 149 in M K and My, respectively, are indicated at a depth of about 
xx20. An excursion 159 in M-m is indicated at a depth of about xx40. It should be 
understood that other formats can be used to present the basic LWD natural gamma ray 
log data. 

5 Fig. 15 shows an example 150 of an azimuthal natural gamma ray LWD log 

presentation. Elemental concentrations of Mr (%), My (ppm) and M-m (ppm) are 
designated with solid, long dashed and short dashed curves 164, 162, and 160, 
respectively, as shown in field 151. Corresponding scales for these concentrations are 
tabulated in field 150 and are typically in % per chart division for M K and ppm per chart 

10 division for My and M-n,. Concentrations M K , My and M-m, obtained from spectra 
summed over a depth interval Ad shown in the field 153, are shown as a function of 
azimuth sector 152, in degrees, for that depth interval. As an example, excursions 156 
and 154 in M K and My, respectively, are shown at an azimuth sector of about 180 degrees 
over the depth interval xx20. As another example, an excursion 158 in Mxh is shown at 

15 an azimuth sector of about 225 degrees over the depth interval xx40. It should be 
understood that other formats can be used to present the basic LWD azimuthal natural 
gamma ray log data. 

While the foregoing disclosure is directed toward the preferred embodiments of 
20 the invention, the scope of the invention is defined by the claims, which follow. 

What is claimed is: 
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